Mean effects of turbulence on elliptic instability in fluids 
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Elliptic instability in fluids is discussed in the context of the Lagrangian-averaged Navier-Stokes- 
alpha (LANS— a) turbulence model. This model preserves the Craik-Criminale (CC) family of 
solutions consisting of a columnar eddy and a Kelvin wave. The LANS— a model is shown to 
preserve the elliptic instability for the inviscid case. However, the model shifts the critical stability 
angle. This shift increases (resp. decreases) the maximum growth rate for long (resp. short) waves. 
It also introduces a band of stable CC solutions for short waves. 

PACS numbers: 47.20.Cq, 47.20.Ft, 47.27.Eq, 47.32.Cc 



Elliptic, or tilting, instability is a fundamental phe- 
nomenon in fluids that results from parametric reso- 
nance. This is the mechanism by which vortex stretch- 
ing creates three dimensional instabilities in swirling two 
dimensional flows. Specifically, the energy in an ellipti- 
cal columnar eddy may be transferred to a propagating 
Kelvin wave Q by this mechanism. A breakthrough in 
the study of this problem occurred when Craik & Crimi- 
nale showed that superposing a columnar eddy and a 
Kelvin wave yields an exact solution of the Euler equa- 
tions. Thus, the elliptical instability can be treated as 
a modulation of the Craik-Crimale (CC) family of solu- 
tions by using Floquet analysis as was first done by Bayly 
The history of discovery and rediscovery of elliptic 
instability in laminar fiuids is reviewed by Kerswell [Q. 
Here we address the mean effects of turbulence on the el- 
liptic instability. The CC family of superposed solutions 
is preserved by the closure model we shall consider. For 
this model, turbulence is shown to enhance the growth 
rates of elliptic instability for Kelvin wavelengths that 
are longer than the turbulence correlation length. Con- 
versely, turbulence is found to suppress elliptic instability 
at the shorter wavelengths and to create a band of stable 
CC fiows with nonzero eccentricities. 

The turbulence model we shall consider is the 
Lagrangian-averaged Navier-Stokes-alpha (LANS— a) 
model |q], whose equations are: 



1^ + (u • V)v + {Vuf ■ V 
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together with V • u = 0. In this model, the mean 
fluid velocity u is related to the mean momentum v via 
the Helmholtz operator (1 — a^A) as v = u — a^Au. 
This Helmholtz filtering of the fluid velocity introduces 
the length scale a as a parameter in the model. The 
LANS— a model preserves the fundamental transport 
theorems for circulation and vorticity dynamics of the NS 
equations. Direct numerical simulations of the LANS— a 
model for forced homogeneous turbulence show it to be 



considerably less computationally intensive than the ex- 
act NS equations while preserving essentially the same 
behavior as NS at length scales larger than alpha |^ . The 
unforced, inviscid Euler— a form of these equations first 
appeared in the context of averaged fiuid models |0] . The 
basic properties of the LANS— a model, its comparison 
with experiment, and its early development are reviewed 
in Ref. ^. See also Refs. ^ for additional results for this 
model. As discussed in Ref. [Q the LANS— a turbulence 
equations formally coincide in the inviscid limit with a 
classic rheological model known as the 2nd-grade fluid 
|]lO| . Thus, the present results for the inviscid elliptic 
instability apply to both the LANS— a turbulence model 
and to rheology of 2nd-grade fluids. 

We construct an exact solution to Eq. (||) with zero 
divergence of the form Uo(x, t) = S{t)x + JJ{t), where 5x 
is the action of the matrix S on the vector x = [x, y, z]'^ 
from the left. The matrix 5 is a time dependent matrix 
with zero trace such that dS/dt+S^ — M{t), where A^(t) 
is a symmetric matrix which contains the contributions 
of U(<). The corresponding pressure po is obtained from 
M{t); see Ref. || for details. We nondimensionalize the 
system using the variables x' = x/l, t' = wt, u' = u/w/, 
v' = v/wZ, a' — a/^/l, where I is a typical length scale 
and u) = i|V x uo|. The resulting equation with the 
prime notation suppressed is Eq. (^ with v replaced by 
v/uj. We construct a second solution to Eq. (|l|) of the 
form Uq + Ui with corresponding pressure po + Pi 7 where 



Ui = ^a(t) sin(/3?/'(x, i)), 
Pi = /^j5ii(t)cos(/37/;(x, t)) 

+ fi^Pi2{t) cos2(/3V(x, t)), 



(2) 
(3) 



t) = k(t) ■x + g{t), and /i and /3 are scaling factors so 
that we can choose the initial conditions |a(0)| = 1 and 
|k(0)| = 1. The unknown phase ip{x,t) and the ampli- 
tudes a(^), pii(t), and puit) are to be determined. Such 
flows which are the sum of a 'base flow' Uq and a 'distur- 
bance' Ui are called a-CC flows. The incompressibility 
condition gives 



ak = 0. 



(4) 
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The evolution equations for the ampHtudes and phase are becomes 



{dt + Sx ■ V)^ + U • k = 0, 
{dt + 5x • V)((l + r)a) + rS'^a + Sa 

-{(3pi, - p^a^a ■ Sk)k = -ii±±^|k|2a, 
pi2-r|a|2 = 0. 



(5) 



(1 + r)i//32 ,. ,o (6) 



(7) 



Here F = a^/3^|kp. Note that the ampHtude scahng /i 
is immaterial. The parameter a couples various terms 
throughout the system. Moreover, this coupling in a ap- 
pears only in the combination F, which is proportional 
to wavenumber-squared. Consequently, this coupling af- 
fects the high wavenumber behavior of the solution for 
a 7^ 0. Equation (||) states that the phase is advected 
with the base flow. Only two free parameters remain 
(F and E^) upon introducing the vorticity based Ekman 
number = v0^ /lo. Without loss of generality, we set 
dg / dt + 'k ■ \J = 0. Denoting the material derivative as 
d* = {dt + Sx ■ V) and taking the gradient of Eq. (||) 
reduces Eqs. (||)-(|6|) to a system of ordinary differential 
equations: 



d((l + r)a) 
dt 



+ (1 + F)5^a 

-2to X a - Pk = -(1 + r)E^\kfa, 



(8) 
(9) 



where P is the coefficient of k in Eq. (||), = x Uq 
is the (normalized) vorticity of the base flow and (5 — 
X a. We eliminate the pressure term by 
taking the dot product of Eq. (^) with k and by using 
da/dt • k — a ■ dk/dt = 5a • k, the first of which follows 
from Eq. ^ and the second from Eq. (^: 

P= ^{(H-F)(5-f5^)a-k + 2ro X a-k}. (10) 

In summary, we have obtained a new exact incompress- 
ible solution to Eq. (|l|). The variables are amplitude 
a{t) and wave vector k(t). Once these are determined, 
the pressure terms follow from Eqs. (|^) and (0). Note 
that uo and uq -I- ui are exact solutions to the nonlinear 
equations, but ui by itself is only a solution to Eq. (|l|) 
linearized about uo. The construction described above 
also can be applied to Eq. (|l|) expressed in a rotating co- 
ordinate system in which an a-CC flow can still be found. 
The effects of rotation will be discussed elsewhere. Fi- 
nally, we emphasize that the operator d/dt -I- S*"^ acting 
on a vector represents the complete time derivative of 
that quantity in a Galilean frame moving with Uq. 

Insight into the dynamics of the problem can be gained 
by examining Eq. (|^) in the asymptotic regimes F <C 1 
and F > 1, where F = a^/g^jkp. (We assume that |k(t)| 
remains bounded and never vanishes.) For F < 1, Eq. (|) 



da rp , 25a • k 

— -f 5^ a = -S^ kra- 2ro x a+ .o k 
at kr 



F ( 2ti7 X a + {(5k • k)a - (nj x a • k)k} 
kM 



+ 0(f2). (11) 

Combined with Eq. ^ , this equation preserves a • k = 
at each order. The 0(1) term in the above equation is 
exactly the expression for the amplitude of the modulated 
traveling wave in the CC flow for the classical Euler and 
NS equations. This, of course, is expected since Eq. (|l|) 
reduces to the NS equations for a = 0. Thus, to leading 
order, the amplitude decays with viscosity, stretches with 
the base shear and rigidly rotates with the vorticity of the 
base flow. For F 1, Eq. (||) becomes 

da „r „ 25a • k, 

— +5^a=-i?^|k|2a+^j^k 

-f ^{(5k.k)a-(tnxa.k)k} + 0(^ij. (12) 

Again, this equation preserves a • k = at each order 
Thus, as F — !■ cx) (corresponding to either a — > oo or 
/3 — > oo), the amplitude no longer rigidly rotates with 
the vorticity of the base flow. 

As an example, we examine the stability of a rotating 
column of fluid with elliptic streamlines whose foci lie on 
the y-axis and vorticity vo — we^ : 



Uo 



wix, L 



-1 + 7 0^ 
1+7 

0; 



(13) 



Here, < 7 < 1 is the eccentricity of the ellipses, and 
the pressure is po — ^^^(1 ~ 7^)(x^ + y"^). Equation (||) 
with 5 = L is analytically solvable: 



k = [sin0cos(i\/l — 7^), 



in(t\/l- 7^), cos 61]^ (14) 



K sm y sm 

where — {1 ~ 7)/(l + 7) and 9 is the polar angle that 
k makes with the axis of rotation. In summary, we have 
a four parameter problem in F, E^, 7, and 9. Eq. (|^) 
has the form da/dt = Af{t)a, where the elements of the 
matrix Af{t) are periodic with period r = 2t^I yj\ — 7^. 
Therefore, the system can be analyzed numerically using 
Floquet theory [ pT| . We compute the monodromy matrix 
P, that is, the fundamental solution matrix with identity 
initial condition evaluated at t = r. Equation (|p) will 
have exponentially growing solutions if max^ |3fi(/9i)| > 
1, where pi^i — 1,2,3 are the eigenvalues of 7^, with 
corresponding Lyapunov-like growth rates given by 



(T = - ln{max |3fi(pj)|}. 

T i 



(15) 
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Thus, we can simulate numerically the solution to Eq. (||) 
over one period and indisputably determine the exponen- 
tial growth rates. We can be certain that at least one of 
the eigenvalues will always be unity because of the in- 
compressibility condition (^) and that the remaining two 
eigenvalues appear as complex conjugates on the unit cir- 
cle or as real valued reciprocals of each other. 

The present investigation considers the case of inviscid 
flow, i.e. = 0. Viscosity, which only slightly modifies 
the inviscid results, will be discussed elsewhere. For flows 
with circular streamlines (7 = 0), the monodromy matrix 
can be analytically computed. It follows from Eq. jl^ ) 
that |k(t)| = 1. Then, T is constant in time (denoted by 
To = o?0^) and Eq. (H) has three linearly independent 
solutions: 



ai(t) = cos(e(t) + </')kii 
a2(t) = sin(e(t)-H0)kii- 
a3(0 = e^. 



- sin(C(t) + 0)k_L2 

COs(C(t) ^ <?!))k_L2 



(16) 
(17) 
(18) 



where i(t) = 2icos6l/(l -I- Eq), k_L2 = [sin t, - cost, 0]^ 
and kj^i = [cos^cost, cosflsini, — sin^?]-^ are vectors or- 
thogonal to k, and is an arbitrary phase. Clearly the 
first two solutions ai and a2 satisfy Eq. (Q). The mon- 
odromy matrix can be constructed from these three so- 
lutions: 

/ cos(^(27r)) cos6'sin(^(27r)) 0^ 

V= \ -sin(C(27r))/cos6' cos(C(27r)) 
\tan6'(l - cos(^(27r))) - sin6'sin(^(27r)) 1^ 

The three eigenvalues are pi^2 = exp(±i^(27r)), = 1. 
All of the eigenvalues lie on the unit circle, from which it 
follows that all solutions in the inviscid case for 7 = are 
stable. The values of cos0 for which \pi\ = l,i = 1,2,3 
are called 'critically stable' and are given by ^(27r) = htt, 
n = 0, ±1, ±2, . . . , corresponding to cos0 = n(l -I- ro)/4. 
At these parameter values an exponentially growing solu- 
tion can appear (together with an exponentially decaying 
one) as 7 increases from zero. Since Fq > 0, the only val- 
ues of interest are n = 0,±1,±2,±3, and, for the case 
a — 0, n — ±4. Bayly |^ argues that the evenness of Fk 
as a function of k implies that the eigenvalues, if real and 
unequal, must be positive. This dismisses the n = ±1 
and n = ±3 choices. The cases n = and n — ±4 pre- 
serve the two-dimensional structure of the base flow and 
thus should be stable under small perturbations in the ec- 
centricity. The remaining value, cos^ = ^(1 + Fq) is the 
critical parameter value at which a(t) suffers exponential 
growth as 7 increases from zero. We conclude that in- 
troducing a preserves the existence of elliptic instability, 
though it shifts the angles at which elliptic instability 
arises to cos 9 = (1 -|- Fo)/2. In addition, for Fq > 1, the 
LANS— a model stabilizes Bayly's elliptic instability. 

Additional understanding of this result emerges by fol- 
lowing the analysis of Waleffe [Q and Kerswell [Q . By 



taking the dot product of Eq. 
all 7) 

d/'i|a|2^ 

= -270102 -f- 



with a, we obtain (for 



dt 



471? ,2 



(19) 



One can determine an exponential growth rate to leading 
order in 7 by inserting the zeroth order solutions for k 
and ai into the right hand side of this equation: 



1 d(i|a|^) 



7 



dt 



[(1 



i0)2sin(2(e+ + 0)) 



(l + cos6')2sin(2(C_ -f 0)) 

2.^..„...M , 27F0 



2(1 - cos^ 0) sin(2t)] 



1 + Fo 



sin(2t), (20) 



where ^± = ^(t) ± t. Upon averaging over a period of 
ai, this quantity will vanish except when ^± = 0, cor- 
responding to cos9 = =pFo/2. The maximum values for 
a will occur at (j) = T7r/4 for ^± = 0, respectively, with 
growth rate 



(3 + Fq)- 
16 



7 + 0(7'), 



(21) 



valid for Fq < 1. Thus, we see that the angle of critical 
stability is again cos 61 = ^(l-|-Fo)/2. Furthermore, we 
see that the maximum growth rate increases as a function 
of Fq due to the Fq dependence of the critical stability 
point up to a maximum of ct = 7 at Fq = 1, after which 
a set of stable solutions emerges in a band of nonzero 
eccentricities. See Fig. |l|. 
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FIG. 1: The growth rate maximized over cos 6 for Eu, — and 
several values of Fo = oP'ff' . The sohd lines are, from bottom 
to top, Fo = 0, 0.1j_0.25, 0.5, 1. The maximum growth rate is 
bounded by Eq. (pl|). The dashed lines, from top to bottom, 
are Fo = 1.25, 2.5, 5.0, 12.5. Notice that for Fo > 1, a stable 
band of nonzero eccentricities appears. 



For nonzero values of 7, we must investigate the system 
numerically. Figure shows the evolution of the critical 
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instability surface as a function of For a^/?^ << 1, 

there is little change in the critical instability surface as 
predicted by Eq. (|ll|). For > 0, all angles of inci- 

dence for the traveling wave are unstable in a neighbor- 
hood of 7 = 1. The maximum growth rate in the j-cosO 
plane increases as a function of a^/?^ and shifts to the 
corner 7=1, cos 9 = 1 by = 0.1. When 

exceeds unity, the flow stabilizes. For a given set of pa- 
rameters (7,cos0), one of following three situations will 
occur as shown in Fig. ^j: the flow is stable for all 
the flow is unstable for < a^(3^ < a\0l and stable for 
0^(3'^ > ajPf; or the flow is stable for < a^/?^ < a^J, 
unstable for 0^(32 < < oi\0l^ and stable again for 




FIG. 2: Surface of a = 0.01 for = 0. The horizontal plane 
is the 7-cosS plane and the vertical axis is ci^ 0^ . Figure 
(a) shows the neutral surface for < 0^0^ < 1 and is an 
expansion of the boxed region in figure (b). For a — 0, the 
critical stability point occurs a.t 9 — n/3, which agrees with 
the classical results. The critical stability point shifts towards 
cos 9 = 1 as a^P^ increases according to cos9 — (1 + q^/3^)/2. 
As a^j3^ exceeds unity, a stable band of rotating flows with 
nonzero eccentricities appears. 

Thus, the LANS— a turbulence model enhances the 



growth rates of the elliptic instability for long waves with 
< 1 while it shifts the angle of critical stability 
along the cusp rising diagonally in Fig. ^. It also stabi- 
lizes the elliptic instability for short waves with a^/?^ > 1 
as seen m Fig. |b. Finally, for any a /^^ ^ 0, this tur- 
bulence model modifies the region in (7, cos 9) parameter 
space where the elliptic instability occurs, as also shown 
in Fig. |. 

In principle, one can now examine the stability prop- 
erties of this new family of exact a—CC solutions of the 
LANS— a model. This would be a secondary stability 
analysis of the rotating base flow. Work of this type 
was carried out by Lifschitz and collaborators for 
the classical CC flows under high-frequency, short wave- 
length perturbations. A similar perturbation analysis for 
the Q-CC flow will be carried out elsewhere. 

The authors are indebted to A. Lifschitz-Lipton for 
stimulating our original interest in CC solutions and 
to the Center for Scientific Computing at Southern 
Methodist University for use of their computing re- 
sources. Furthermore, BF thanks the Theoretical Di- 
vision at the Los Alamos National Laboratory for their 
hospitality. The solutions to a(t) were simulated using 
the variable coefficient ODE solver DVODE 
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